Geometric Scaling from DGLAP evolution 
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We show that the geometric scaling of the total virtual photon-proton cross section data can be 
explained using standard linear DGLAP perturbative evolution with generic boundary conditions 
in a wide kinematic region. This allows us to single out the region where geometric scaling may 
provide evidence for parton saturation. 
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The observation of geometric scaling [1] in ep deep in- 
elastic scattering (DIS) data has attracted considerable 
interest because it is widely interpreted as evidence for 
parton recombination and saturation Q ■ Geometric scal- 
ing (GS) is the statement that the total "f*p cross section 
aJ oi f, which is a priori function of two independent vari- 
ables — the photon virtuality Q 2 and the Bjorken vari- 
able x — only depends on the variable r = Q 2 /Q 2 (x), 
where the so-called saturation scale Q 2 (x) depends non- 
trivially on x, with dimensions given by a fixed reference 
scale Qq. 

The presence of recombination effects in the HERA 
data would have dramatic effects, because these data 
dominate the determination of parton distributions, 
which are necessary for the computation of LHC pro- 
cesses Q. Available parton fits do not include these ef- 
fects and would thus fail to provide reliable predictions 
at the LHC. It is the purpose of this paper to ascertain 
whether this is actually the case, and more in general, 
in which kinematic region GS may or may not provide 
evidence for saturation. 

Evidence for GS is provided by the scaling plot of the 
reduced cross section a] f vs. the scaling variable r: 
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FIG. 1: Data sample. Lines of constant r are shown, with 
lnr = t — A£ (dashed) and lnr = t — Ay^ (solid) and Q\ = 
1 GeV 2 . 
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where t = \hQ 2 /Qq, and the saturation scale t s = 
In ®qi^ was originally 1] chosen as t s = A£, and more 

recently @, H also as t s = A^f, with £ = ln(l/x). To 
test whether this behaviour is compatible with standard 
DGLAP perturbative evolution, in Figs. EE] we show a 
scaling plot of a^ f computed using the double asymp- 
totic scaling (DAS) approximation to leading-order (LO) 
DGLAP evolution [5] of a constant boundary condition. 
Namely, we take 



FIG. 2: Geometric scaling for the data of Fig. 1 with lnr = 
t - A£ and Q 2 , = 1 GeV 2 . Only data with Q 2 > 1 GeV 2 are 
in the DAS and improved DAS curves. The DAS curves are 
offset for clarity. 



The DAS approximation to LO DGLAP evolution is 
quite accurate in a wide kinematic region Q 2 ^ 10 GeV 2 , 
x ^ 0.1, based on approximating the LO anomalous di- 
mension as 
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We take it as representative of a situation where satura- 
tion is certainly absent. The GS properties of the cross 
section Eq. ([2]) are compared to those of the data by 
plotting both vs. r, with the best-fit phenomenological 
values of Ref. Q A = 0.321 (Fig.© or A = 1.621 (Fig.©, 
over the grid of (x, Q 2 ) values shown in Fig.[TJ The data 
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FIG. 3: Same as Fig. d but with Inr = t - Ap- 



points arc obtained from a very accurate neural network 
interpolation to world DIS data 0]. The kinematical re- 
gion is that for which experimental data are available. 

Figs.EElshow that for r > 1 (Fig. 2) or r > 0.1 (Fig. 3) 
the GS properties of the DAS solution are almost as good 
as those of the data, and become as good or better with 
a minor improvement, to be discussed below. It follows 
that saturation is by no means necessary for geometric 
scaling. This may seem surprizing given that the DAS 
solution Eq. @ appears to violate GS. However, we now 
show that approximate GS is in fact a general property 
of solutions to the DGLAP equation. It is already known 
that if GS is imposed as a boundary condition at some 
low scale, it is preserved by both BFKL [2] or DGLAP Q 
linear evolution to higher scales. Here, we show instead 
that GS is generated by linear DGLAP evolution itself, 
irrespective of the choice of boundary condition. 

Consider first the fixed-coupling case. The general 
DGLAP solution for any anomalous dimension 7 is 



— G (N)exp[N^ + 1 (a s ,N)t} 7 (5) 
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where Gq(N) is a suitable boundary condition. For large 
enough £, the integral can be evaluated in the saddle 
point approximation. The saddle condition is 



The cross section becomes 
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up to terms which are not enhanced as £ — > 00. 

Geometric scaling follows expanding t about the satu- 
ration scale Eq. (JTJ) t s = A£: 
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If we choose a value of A such that /(A) = the cross 
section Eq. ([5]) manifestly displays geometric scaling. It 
is apparent from Eq. Q that this value exists if 7 is 



the DGLAP anomalous dimension, either at fixed per- 
turbative order or resummed at small x using the BFKL 
formalism, or indeed for any reasonable shape of 7. 

This argument is in fact quite close to that of Ref. [|J , 
due to the fact that the DGLAP solution can equivalcntly 
be written in "dual" form Q as 

f c+tQC dM - 

G(t, = / — G (M) exp [Mt + x (a s , M)£] , 

(9) 

where the kernel \ is related to 7 by 

X [a s , 1 {a s ,N)]=N, (10) 

and Go (M ) is determined in terms of the boundary con- 
dition Gq(N) and the anomalous dimension 7. Evaluat- 
ing the integral by saddle point and then Taylor ex- 
panding reproduces the argument for geometric scaling 
of Ref. [2]. However, Eq. ([8]) shows that the "saturation" 
assumption of Refs. 0, 0] that the boundary condition 
satisfies GS is redundant: rather, GS follows from the 
existence of A such that /(A) =0 in Eq. ©. This is a 
generic property of perturbative evolution. Equation © 
can be equivalently viewed as the solution to the BFKL 
or DGLAP equations, and our conclusion applies to both. 

We conclude that GS holds for the solution to the 
DGLAP equation at the fixed coupling level, which ex- 
plains the GS properties of the DAS solution Eq. @, 
Fig. [2j this solution is derived with running coupling, but 
in practice (see Fig.[T]), the value of t along fixed r curves 
is almost constant in the data region. It follows that 
■i^ ln[(i + to) /to] ~ o-s(Qo) t) which in turn implies that 



Eq. ^ holds with a ~ y/£Poa s t and p w '£/ '((3 a s t), 
which coincides with the result found using Eq. ^ in 
the approximation Eq. |7[). Higher order terms in this 
expansion lead to GS violations, proportional to powers 
of a s (Q^)t. The combined effect of GS violations will be 
discussed in Fig. 4 below. 

A running coupling form of GS can also be derived Q 
directly for the cross section Eq. ([2][3]). At the running 
coupling level, we can neglect the variation of hit /to in 
a and p in comparison to the scale dependence of Q~ 2 
in Eq. |[2}. Then, the DAS solution {3J only depends on 
a/J, and the cross section $2$, consistently neglecting the 
variation of In £ in comparison to the variation of £, is a 
function of the scaling variable t — A-\/£- 

Note that, unlike the fixed coupling GS in terms of 
t — A£, Eq. ((8]), which holds for a generic anomalous 
dimension 7, this running coupling GS depends on the 
particular form of the anomalous dimension Eq. Q , and 
specifically on the fact that it has a simple pole at N = 0. 
However, this running coupling GS can also be obtained 
using the running-coupling version of Eq. ([9]) l3, 11 1 
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FIG. 4: Values for A determined from Eq. (fT2)l . 
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FIG. 5: The quality factor Q computed for the points of 
Fig. [2] with Q 2 > 25 GeV 2 . The solid curve is in each case 
the result of a gaussian fit. 



which holds whenever the kernel x Eq. (|T0|) is linear in 
a s . This is the case if we only retain the simple pole in 
the anomalous dimension 7 Eq. (j4|) , but also for a generic 
leading-order BFKL kernel, as discussed in Ref. [12J. The 
argument leading to GS Eq. (JSJ) can now be repeated: the 
only difference is that the saddle point depends on t j 
which leads to the form t s = A-^/f of the saturation scale. 

The arguments so far involved several approximations. 
To assess their accuracy, we express the cross section as a 
function of In r = t — A£ and the orthogonal combination 
C = t + A£. Geometric scaling is the statement that a] Q f 
is independent of £: 
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By letting A depend on £ and i, we can view the condi- 
tion (fT2"J) as an implicit equation for A(£,£). Geometric 
scaling holds to good approximation if the solution A(£, t) 
to Eq. (|12p is approximately constant in £ and t in the 
kinematic region of interest. We see from Fig. 0] that 
the value of A is almost constant everywhere, except at 
low Q 2 25 GeV 2 . At low Q 2 , large x there are no data 
(see Fig. [1]); while the low Q 2 , low x region, where the 
DAS solution is not applicable, shall be discussed below. 
We would now like to test whether instead in the region 
Q 2 <; 25 GeV 2 the approximate GS displayed in Fig. 2] is 
sufficient to explain the GS of the data. To this purpose, 



we use the "quality factor" Q(X) which was introduced in 
Ref. [3| as a measure of the scaling quality. The optimal 
value of A is that which maximizes Q, and GS is better 
if Q(X) is larger. In practice, the optimal value of A and 
the uncertainty on it are determined by fitting a gaus- 
sian form to Q(X). The values used to produce Figs.OO 
taken from Ref. [J], were determined thus. 

In Fig. [5] we display the quality factor Q(X), computed 
for all points with Q 2 > 25 GeV 2 , both for the data and 
the DAS solution, as well as the result of a gaussian fit. 
We see that GS is actually rather better for the DAS 
solution than for the data. However, we also see that 
the optimal value of A for the DAS solution is somewhat 
larger. This explains why the GS properties of the DAS 
solution in Fig. [5] actually look a bit worse: the optimal 
value of A is not quite the same for the data and for the 
DAS solution. 

To explain this difference, we note that for 

medium-large Q 2 the DAS approximation can be sub- 
stantially improved by including the contribution of the 
smaller eigenvector of the anomalous dimension matrix, 
whereby in Eq. © F 2 = (j/p)G + G, where G(£,t) = 
kexp(—5a/p), with S = 16n//(27/?o) and k — 0.16 is de- 
termined from a fit to the data. The GS plot for this 
improved DAS solution is also shown in Fig. [2l and the 
corresponding quality factor is displayed in Fig. [5j GS 
deteriorates slightly for Q 2 > 25 GeV 2 , but remarkably 
it now holds for all data of Fig. [2l Also, the optimal value 
of A extracted from the data and the improved DAS so- 
lution now agree. This means that linear leading-order 
perturbative evolution, as embodied by the (improved) 
DAS solution, can actually predict the optimal choice of 
saturation scale Eq. {]]). Indeed, a gaussian fit to the 
quality factor for the improved DAS solution on all points 
of Fig.[2](Q 2 > 10 GeV 2 ) gives A = 0.32±0.05, in perfect 
agreement with the value A = 0.32 ± 0.06 determined in 
Ref. [4| from the data. We take this as very strong ev- 
idence that GS for Q 2 Jj 10 GeV 2 follows from purely 
linear perturbative arguments. 

A similar analysis based on Eq. (fT2"|) , but with the run- 
ning coupling form of In r = t — A-y/C an d C — t + -^VC 
leads to the same conclusion. In particular, in this case 
we predict A = 1.66 ± 0.34 to be compared to the exper- 
imental value A = 1.62 ± 0.25 of Ref. 

It remains to be understood why the data still display 
GS even at low Q 2 where the DAS solution becomes un- 
realiable. Figure 2] suggests that the effective value of A 
which characterizes perturbative evolution starts grow- 
ing significantly for Q 2 ^ 10 GeV 2 . Because GS is never- 
theless seen in this data region (see Fig. [2][3|), one might 
conclude that there is some evidence for saturation there. 

However, so far we have only used pure leading-order 
DGLAP evolution, which fails in this region because it 
does not resum small x logarithms. Before conclud- 
ing that a saturation-based approach is necessary, we 
should address the issue of small x resummation in the 
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framework of linear perturbative evolution. The small 
x resummation of DGLAP evolution has been recently 
performed, based on a suitable matching of the BFKL 
and DGLAP solutions 14, [lij]. For our present pur- 



pose, it is enough to consider the asymptotic small x be- 
haviour of these matched solutions, which is essentially 
determined ll|, [14J by a quadratic approximation to a 



running coupling BFKL evolution kernel. 

As is well known , this leads to a solution written in 
terms of Airy functions Ga(N, t) if the kernel is linear in 
a s (Bateman functions [14j if the nonlinear dependence 
is retained). One can then extract ll| an anomalous 
dimension 



j A (a s (t),N) = -\nG A (N,t). 



(13) 



The asymptotic small x behaviour of DGLAP evolution 
at the resummed level is controlled by the rightmost sin- 
gularity of 7A(a s (i), N) Eq. (fT3|) . This singularity turns 
out to be a simple pole, located at N = N (t). 

Neglecting the weak [ll[ scale dependence of N a the 
predicted asymptotic small x behaviour at the resummed 
level is 
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This behaviour should hold in a region where Q 2 is large 
enough for some resummed linear perurbative evolution 
from a low-scale boundary condition to have taken place, 
say 5 < Q 2 < 10 GeV 2 . In this region, we thus get GS 
with A = No. Typical values of A*o from resummed linear 
perturbative evolution are 0.1 < N Q < 0.3 Q. The GS 



properties of the unresummed DAS solution are thereby 
extended down to scales of order 5 ^ Q 2 ^ 10 GeV 2 . 

The scale dependence of N can be kept into account 
by determining it in an expansion in powers of a 2 / 3 (t) (llj 



N (t) =ca s (t) 



where c = x 
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and k = 



2^7 mp X q {a s ,M)\ M =M a 
parametrize the quadratic BFKL kernel x q ( a s, M) at its 
minimum M = M . 

Substituting this in Eq. (frjj and expanding we see that 
the asymptotic form of the cross section is constant on 
the curve 
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This corresponds to the previously discussed running 

coupling form of the saturation scale, with A = \J^- 

Realistic values 1 < c < 2 give 1.2 < A < 1.7. This 
implies GS along this saturation line, and approximate 
GS in the proximity of it. Of course, if we go very far 



from this region we end up in the large Q 2 region which 
we have already discussed. 

Hence, thanks to small x resummation the growth 
Fig. 0] of A at small x for 5 £ Q 2 ^ 10 GeV 2 is replaced 
by the (almost) constant value Nq, thus extending GS to 
this region, with t s — A £ or t s = ^/47rc//3o\/£- 

In conclusion, we have shown that for Q 2 <; 10 GeV 2 
standard linear leading-order DGLAP perturbative evo- 
lution explains geometric scaling, and in fact predicts 
the value of the constant A which characterizes the "sat- 
uration" scale Eq. ([I}. Small x resummation of the 
linear evolution equation extends the region where GS 
is expected to values of Q 2 which are lower, but still 
within the perturbative region. For yet lower values of 
Q 2 5 GeV 2 , geometric scaling, which is observed in the 
data, cannot be explained using linear perturbation the- 
ory. This is the region r 5s 0.1 (fixed-coupling t s , Fig. ^ 
or t ^ 0.01 (running-coupling t s , Fig. [3]), where the GS 
plot flattens out. In this region, geometric scaling may 
provide genuine evidence for parton saturation. 
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